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ABSTRACT 

We present two new in situ core accretion simulations of Saturn with planet formation timescales 
of 3.37 Myr (model SO) and 3.48 Myr (model SI), consistent with observed protostellar disk lifetimes. 
In model SO, we assume rapid grain settling reduces opacity due to grains from full interstellar val- 
ues (Podolak 2003). In model SI, we do not invoke grain settling, instead assigning full interstellar 
opacities to grains in the envelope. Surprisingly, the two models produce nearly identical formation 
timescales and core/atmosphere mass ratios. We therefore observe a new manifestation of core accre- 
tion theory: at large heliocentric distances, the solid core growth rate (limited by Keplerian orbital 
velocity) controls the planet formation timescale. We argue that this paradigm should apply to Uranus 
and Neptune as well. 

Subject headings: planets and satellites: formation — planets and satellites: individual (Saturn) 



1. INTRODUCTION 

The discrepancy between the observed lifetimes of pro- 
tostellar disks (2-3 Myr; Haisch et al. 2001) and the 
length of time required for planet formation by core ac- 
cretion (> 8 Myr; Pollack et al. 1996) has long presented 
a problem for planet formation theory. However, in 2005, 
two new models of planet formation showed that the 
core accretion-gas capture process could form gas giants 
within 2.5 Myr. Hubickyj et al. (2005) modeled Jupiter's 
formation in situ at 5.2 AU and found that Jupiter could 
grow from a 0.1M® core to its present mass in 2.2 Myr. 
The core accretion models of Alibert et al. (2005) allowed 
both Jupiter and Saturn to form concurrently within 
2.5 Myr in disk of mass 0.035 M < M disk < 0.05 M . 

Each model employed a different approach in order to 
speed up giant planet formation. Alibert et al. added 
Type I migration to the core accretion model: the in- 
ward motion of the protoplanets allows them to receive 
a fresh supply of planetesimals and gas as they move 
into undepleted regions of the solar nebula. The accre- 
tion rate is then no longer limited by the rate at which 
the protoplanet's Hill sphere expands, so the giant planet 
formation timescale decreases by up to an order of mag- 
nitude. These models require proto-Jupiter to have an 
initial semimajor axis of a > 9.2 AU and proto-Saturn to 
begin forming at 11.9 AU, and migrate to their current 
positions on a ~ 2 Myr timesc ale. 

One notable feature of the lAlibert et al.l (|2005h core 
accretion model is that it successfully predicts the heavy 
metal content of Jupiter and Saturn's atmospheres ac- 
cording the clathrate hydrate trapping theory of Luninc 
& Stevenson (1985). However, the Type I migration rate 
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is a free pa rameter: Alib e rt et al. note that the analyti- 
cal work of lTanaka et al.l (|2002f l predicts migration rates 
far too large to be consistent with the observed frequency 
of extrasolar plane ts. The authors get a round this prob- 
lem by scaling the iTanaka et al.l (|2002l ) by an arbitrary 
factor fx, where < fx < 0.03. Finally, there is one 
more basic assumption underlying the model, which is 
that the gas/solid ratio beyond the ice line is constant at 
G/S = 70. 

The in situ planet formation models of Hubickyj et 
al. (2005) decrease Jupiter's formation time by requiring 
that grains quickly settle to the bottom of the proto- 
planetary envelope, where they are destroyed by subli- 
mation (Podolak 2003). Assuming grain settling lowers 
envelope opacity due to grains to ~ 2% of the interstel- 
lar value, the gas can contract efficiently and make way 
for new material entering the protoplanet's Hill sphere. 
The most important free parameter in the Hubickyj et 
al. models is the solid surface density of planetesimals in 
the planet's feeding zone: since the protoplanet doesn't 
move through the disk, it requires a feeding zone with 
Ssoiid ^ 10 g cm" 2 in order to attain the ~ 15M core 
necessary for accreting a massive gaseous envelope (Pa- 
paloizou & Nelson 2005). 

The availability of new calculations of solid surface 
density as a function of heliocentric distance in the solar 
nebula (Robinson et al. 2008; hereafter Paper 1) raises 
the possibility of extending in situ core accretion simula- 
tions to include Saturn. Since the Nice model of Tsiganis 
et al. (2005) predicts that proto-Saturn migrated 1 AU 
outward at most, we consider the in situ approximation 
reasonable when applied to Saturn. By providing the- 
oretically and observationally motivated values for solid 
surface density S S oiid, the Paper 1 results move core ac- 
cretion simulations away from parameter studies and to- 
ward determinism. 

Since the solar nebula had to be capable of forming 
both Jupiter and Saturn concurrently (and of course 
Uranus and Neptune, the formation of which we will ex- 
amine in future work), we first assess the ability of the 
Paper 1 solar nebula model to produce Jupiter. In Paper 
1, the value for £ so iid after 10 5 yr of solar nebula evolu- 
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tion is 13.2 g cm -2 . Adopting the relationship between 
solid surface density and Jupiter formation time at 5 AU 
calculated by Robinson et al. (2006), 

kform / ^solid \ \ 

lMyr ~ V 25 -°g cm_2 / ' 
we find that the Paper 1 r esults allow Jup i ter to form 
in 2.5 Myr. Built in to the iRobinson et alj (|2006l ) scal- 
ing relation is the assumption that efficient grain settling 
leads to protoplanetary envelope opacities of ~ 2% those 
of interstellar grains. In this Letter, we will relax this as- 
sumption and also investigate the limiting case of 100% 
interstellar grain opacity with respect to Saturn. 

If the fiducial disk from Paper 1 can form Saturn within 
2-3 Myr, we will have successful core accretion models of 
the two gas giants forming near their present positions in 
a gravitationally stable disk (see Paper 1 for a discussion 
of the solar nebula dynamics). 

In ij^l we describe our theoretical treatment of the core 
accretion process. In Jll we discuss the results of our sim- 
ulations, with special emphasis on formation timescale, 
effect of atmospheric opacity, and core/atmosphere mass 
ratio. We present our conclusions in |J4] 

2. THEORETICAL TREATMENT OF CORE ACCRETION 

We use the theoretical model of planet formation de- 
scribed by Laughlin et al. (2004) to model the core ac- 
cretion and gas capture of proto-Saturn. Initially, a pro- 
toplanetary core of mass M® is embedded at Saturn's 
heliocentric distance, 9.5 AU, in a viscously evolving disk 
of age 1.5 x 10 5 yr, surrounding a T-Tauri star of mass 
IMq. We assume that by 1.5 x 10 5 yr, the available 
dust has formed 100 km planetesimals that are invul- 
nerable to gas drag (Weidenschilling 1977): planetesi- 
mal orbits are modified only by interactions with proto- 
Saturn. Gas temperature and density are regulated by 
viscous evolution of the solar nebula. We use the time- 
evolving temperature and density at 9.5 AU, beginning 
at t = 1.5 X 10 5 yr, as calculated in Paper 1. 

The contraction and buildup of protoplanetary cores 
and their gaseous envelopes embedded in our model 
evolving disk are computed with a Henyey-t ype code 
Heny ey et al. 1964). Following the argument of lPodolaki 
2003) 'that grain settling in the protoplanetary enve- 
lope would reduce envelope opacity where grains exist, 
we ado pt grain opaci ties of 2% of the interstellar values 
used in lPollack et all (|1996h in our fiducial model, which 
we will call SO. However, in order to assess the effect 
of envelope opacity on Saturn's formation timescale, we 
present a second core- accretion simulation, SI, using full 
interstellar grain opacity. 

We use a core accretion rate of the form 

— Jj^L — Ci7rE solid i? c i?/ l f7 (2) 
at 

(Papaloizou & Terquem 1999), where E S oiid is the surface 
density of solid material in the disk, f2 is the orbital fre- 
quency at 9.5 AU, R c is the effective capture radius of the 
protoplanet for solid particles, Rh — a[M P ianot/(3M*)] 1 / 3 
is the tidal radius of the protoplanet (where a is the semi- 
major axis of the protoplanet's orbit), and C\ is a con- 
stant near unity. 

The outer boundary conditions for the protoplanet in- 
clude the decrease with time in the background nebu- 
lar density and temperature. During the late phase of 



planet growth, when planetesimals may be ablated by 
the massive envelope before reaching the core, we con- 
sider a planetesimal captured if it deposits 50% or more 
of its mass in the envelope. At this stage, we invoke 
the sinking approximation and assume the ablated plan- 
etesimal debris sinks rapidly to the planet core without 
leaving remnants in the envelope. 

3. RESULTS 

In both simulations (SO and SI), we start the core ac- 
cretion model with midplane temperature, gas density 
and solid surface density from the solar nebula model of 
Paper 1. This model has two key features favoring planet 
formation that are missing from passive disk models: (1) 
viscous stresses drive the initial E cx R~ 3 / 2 surface den- 
sity profile toward uniformity, so that Saturn's feeding 
zone gains mass during the first 5 x 10 4 yr of disk evo- 
lution, and (2) the presence of hydrated ammonia ice at 
the snow line increases the solid surface density by 7% 
over the standard water ice-rock-refractory CHON mix- 
ture. We use a starting solid surface density of 8.6 g cm~ 2 
(Paper 1) which decreases with time as proto-Saturn cap- 
tures planetesimals. 

Figure la shows the growth of Saturn from a core of 
IMq to its present-day mass of 95M®. Solid lines cor- 
respond to the SO model (in which grains quickly settle 
and sublimate, reducing their contribution to envelope 
opacity by an assumed factor of 50), and dashed lines 
show the SI model (in which grains stay in the envelope 
and opacity due to grains takes on the full interstellar 
value) . An important property of our model, and one of 
two key results of this Letter, is Saturn's formation time: 
the planet attains its current mass in only 3.37 Myr for 
model SO and 3.47 Myr for model SI. These are the first 
in situ core accretion models of Saturn with formation 
times within 1 a of observed protostcllar disk lifetimes. 

The second key result is that Saturn's formation time 
is nearly independent of the assumed grain opacity in the 
envelope. This surprising result occurs because Saturn's 
core growth rate is limited by the Keplerian speed in 
the feeding zone: a slow-moving core takes a long time 
to encounter and capture planetesimals (see Equation 
[2]), and may never reach isolation mass. Whereas during 
Jupiter's formation early core buildup takes only 0.5 Myr 
and the gas contraction phase dominates planet growth 
(Hubickyj et al. 2005), in Saturn's case, core growth lasts 
until the planet reaches its final mass. Only in the last 
6 x 10 5 yr of growth does Saturn possess a gas mass of 
>lM e . 

Given that the planet has a negligible envelope mass 
throughout most of its formation, it stands to reason 
that envelope opacity would not exert much influence on 
Saturn's formation timescale. The idea that solid growth 
rate controls Saturn's formation is consistent with the 
planet's high core/ atmosphere mass ratio: 9 — 23%, as 
opposed to < 3% for Jupiter (Saumon & Guillot 2004). 

3.1. Core/ 'Atmosphere Mass Ratio 

If all accreted planetesimals and their ablated debris 
end up in the solid core, as is assumed in our model, 
Saturn's total core mass reaches 44M® in model SO and 
54M® in model SI. Based on gravitational moment mea- 
surements and internal structure modeling, Saumon &z 
Guillot (2004) place Saturn's core at 9-22M® and total 
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heavy element content at 13-28Mq. Our model Saturn 
has a heavy element mass that is too high by at least 
a factor of 1.6 (model SO). We note, however, that the 
carbon enrichment in Saturn's atmosphere has recently 
been revised upward to (C/H)/(C/H)q = 7 (Flasar et 
al. 2005), which may allow for a larger heavy element 
inventory than previously thought. 

One method of reducing the core mass, though not the 
total heavy element mass, is to account for planetesimal 
disruption not only by ablation of debris (|J2) , but by sub- 
limation of volatiles. Our planetesimals are 50% H2O by 
mass, 6% NH3, and 1% other ices, such as HCN and H2S, 
for a total of 57% ice. These ices could sublimate cither 
during infall through the envelope, analogously to mete- 
orites in Earth's atmosphere, or upon impact with the 
solid core. Helled et al. (2008) found that sublimation of 
volatiles during the collapse of a giant planet formed by 
disk instability (Boss 2005) is efficient: refractory silicate 
grains sediment to form a core, while ices remain in the 
planet atmosphere. 

Tingle et al. (1991) tested the survival of volatiles 
experiencing hypervelocity impacts (v ~ lkms -1 ) by 
shocking samples of the Murchison meteorite with pres- 
sures up to 36 GPa. They found that 70% of organic 
and organosulfuric material, including H2S, sublimates 
upon experiencing an impact with v > 1.5 km s -1 . The 
volatiles in Saturn-building planetesimals are not likely 
to survive a high- velocity impact and subsequent contact 
with the hot protoplanetary core, T > 3000 K, in solid 
form. If we assume all accreted ices undergo a phase 
transition from solid to gas, Saturn's core mass drops to 
19M e (SO) and 23M ffi (SI), which are near the range 
determined by Saumon & Guillot (2004). However, the 
total heavy element/hydrogen mass ratio in the planet 
is still higher than observed: Saturn is 50-70% hydrogen 
by mass, whereas our model predicts a hydrogen mass 
fraction of only 42% (SO) and 31% (SI). 

Another possible way to reduce Saturn's heavy element 
mass without slowing the planet's growth is to cut off 
solid accretion midway through planet formation. This 
approach simulates the effect of another embryo compet- 
ing for planetesimals. Although our solid accretion rates 
are calculated for the monarchic growth paradigm, in 
which planetesimal dynamics are determined by proto- 
Saturn only, this scenario is an approximation: Kokubo 
& Ida (1998) predict an oligarchic planet formation epoch 
with competing embryos spaced ~ 10 Hill radii apart 
that lasts for ~ 1 Myr. 

Hubickyj et al. (2005) tested the effect of a core ac- 
cretion cutoff on Jupiter's formation and found that as 
long as the core has mass > 10M® before solid accretion 
ceases, the planet reaches hydrodynamic gas accretion 
even more quickly than when solid accretion continues 
unchecked: 0.78 Myr vs. 2.22 Myr. Continuous, late- 
stage planetesimal accretion slows planet formation by 
depositing kinetic energy in the protoplanetary envelope 
and inhibiting gas contraction. 

The competing embryo scenario holds promise for 
bringing Saturn's core/ atmosphere mass ratio into agree- 
ment with observations while still retaining the quick 
formation time. However, it is a double-edged sword: 
oligarchs with overlapping zones of gravitational influ- 
ence can increase the RMS planetesimal eccentricity, 



(e 2 ) 1 / 2 , far more efficiently than a single monarch core, 
decrea sing the gravitational focus ing ability of all em - 
bryos. iThommes et~aH (|2003l ) and llda fc Makinol (|1993l ) 
find that significant planetesimal stirring can occur when 
embryos are betwe en 10~ 5 and 10~ 2 M e . Indeed, 
iFortier et al.l (|2007h replaced Eq uation [2] with the an - 
alytical oligarchic growth rate of llda fc Makinol (|1993f ) 
in in situ core accretion simulations of Jupiter and found 
formation times of 10-20 Myr. 

Thommes et al. (2008) created a self-consistent planet 
formation model that included gas disk evolution, planet- 
disk interactions (including gap opening and gas accre- 
tion onto solid cores), and planet-planet interactions. 
Gas giants with Jupiter-like core/ atmosphere mass ra- 
tios consistently emerged in disks with mass M^isk ^ 
0.06M Q . Even a modest amount of migration, as oc- 
curred in the Solar System, appears to compensate for 
the inefficiency of in situ oligarchic growth. The true Sat- 
urn formation scenario probably involved several ~ 1M 8 
embryos spanning 9-12 AU, some outward motion of 
proto-Saturn, and possible competition with Neptune, 
forming near 12 AU (Tsiganis et al. 2005). 

3.2. Isolation Mass 

One new feature of both models, SO and SI, is that Sat- 
urn does not need to accrete all the planetesimals in its 
feeding zone to reach its current mass. Figure lb shows 
the remaining solid surface density, not incorporated in 
Saturn's core, as a function of time. (The upturn near 
the end of the simulation results from the rapid expan- 
sion of the planet's feeding zone during hydrodynamic 
gas accretion and the assumption that planetesimals are 
always uniformly distributed, which is unphysical if the 
planet grows more quickly than planetesimal redistribu- 
tion can occur.) In both models, there are still 2 g enr 2 
of solids left in Saturn's feeding zone at t = 3.4 Myr. This 
means the core never reaches isolation mass. Following 
Lissauer (1993), we calculate 

M iso = 0.0021 E s 3 ^ d a 3 = 58M e , (3) 

whereas Saturn's core attains a mass of 44M® in model 
SO and 54M e in model SI. 

Both Hubickyj et al. (2005) and Pollack et al. (1996) 
found that giant planets pass through a lengthy plateau 
phase in which both solid and gas accretion rates are 
low and planet mass changes very little with time. This 
phase begins when the planet core nears isolation mass, 
and its duration is regulated by gas contraction efficiency. 
Examining Figure la, we see that Saturn never experi- 
ences this mass plateau: both solid and gas accretion 
rates increase with time. 

With extra solid and gas mass remaining in the disk 
after 3.4 Myr, why would Saturn not continue to grow? 
We arbitrarily stop the simulation once Saturn reaches 
95M0, but our disk could quickly form a more mas- 
sive planet. Explaining Saturn's low mass in the con- 
text of the disk models in Paper 1 may require either 
another nearby embryo competing for both solids and 
gas-perhaps Neptune, forming near 12.5 AU according 
to the Nice model (Tsiganis et al. 2005) — or a mech- 
anism for disk dissipation that begins after a few Myr. 
Photoevaporation (Alexander et al. 2006) is one possi- 
bility: the disk near 10 AU begins to be disrupted by 
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Fig. 1. — a (Left): Mass of Saturn as a function of time. Solid lines represent model SO (reduced grain opacity) and dashed lines 
represent model SI (full grain opacity). The black curves show the total planet mass, the blue curves show the solid mass only (presumed 
to be concentrated in the core), and the red curves show the gas mass. Saturn reaches its current mass, 95A/0, in 3.4 Myr. b (Right): 
Remaining solid surface density, not yet accreted by Saturn, as a function of time. Solid line shows model SO and dashed line shows model 
SI. The sharp upturn at 3.3 Myr in SO is due to the rapid expansion of Saturn's feeding zone when hydrodynamic gas accretion begins. 



ionizing radiation at 2-3 Myr. We also require a mech- 
anism to populate the Oort cloud with 40M^ of comets 
(Weissman 1996). Saturn scattering the 4 — 14M e it does 
not accrete into the Oort cloud would be a good start, 
and Jupiter could contribute more material during the 
late stages of its formation and/or migration. 

Is the concept of isolation mass relevant for any giant 
planet except Jupiter? From Equation [3l we see that 
Mj so increases with radius if a < 2 in the power-law sur- 
face density profile £ oc R~ a . Extrapolating from the 
simulation presented here and the near-flat solid surface 
density profile seen in Paper 1, we propose that Saturn, 
Uranus and Neptune do not reach isolation mass. Due 
to their low Keplerian orbital speed, planet cores in the 
outer solar nebula experience a lower planetesimal colli- 
sion rate than proto-Jupiter. Instead of their formation 
timescale being governed by gas contraction efficiency, 
the solid accretion rate, mediated by Q and S, is the crit- 
ical factor. The hypothesis that the three outer planets 
never reach isolation mass is consistent with the plan- 
ets' bulk composition: all three have solid/gas ratios far 
higher than Jupiter. 

4. CONCLUSIONS 

Based on the solid surface densities calculated in Paper 
1, we have created core accretion models of Saturn with 



a formation timescale of 3.37-3.48 Myr, which is com- 
patible with observed protostellar disk lifetimes. Unlike 
previous studies of Jupiter (e.g. Hubickyj et al. [2005]), 
we find that grain opacity in the protoplanetary envelope 
has virtually no effect on Saturn's formation timescale. 
The same solar nebula model that was the basis for 
simulations SO and SI is capable of forming Jupiter in 
2.5 Myr, assuming reduced grain opacities. 

Finally, this is the first in situ core accretion model 
of a gas giant that lacks a plateau phase, in which the 
planet's mass remains nearly constant at approximately 
the isolation mass. We postulate that the low orbital 
speeds far from the sun prevented Saturn, Uranus and 
Neptune from accreting solids efficiently enough to reach 
isolation mass. These planets never experienced a gas- 
only accretion phase, as did Jupiter, and consequently 
have much higher core/ atmosphere mass ratios. 
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